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Abstract 


Collisional  efficiencies  for  vibrational  energy  transfer  have  previously 

These  quantities  have  been 
examined  here  as  a  function  of  molecular  complexity,  critical  reaction 


been  related  to  the  quant i t i es, <AE>,  and  <AE> 

a  a 


threshold,  E^,  temperature  and  energy  transfer  step  size.  Model  calculations 
were  made  for  the  decomposition  of  nitryl  chloride  and  for  the  i somer i zat i ons 
of  methyl  isocyanide,  cyclopropane  and  eye  I oheptatr I ene .  For  large  step  size 
and  low  temperature,  <AE>^ I  I  =  -<AE>^;  at  high  temperatures  up  transitions 


dominate  and  <AE>  =  <AE> 


An  inversion  temperature  T^  is  predicted  when 


up 


down 


all  ■  ' d ' 

;  this  temperature  is  quas i - i ndependent  of  step  size  and  dependent 


on  molecular  energy  and  molecular  complexity,  Tj  is  used  to  compute  an 
elFective  temperature  which  can  be  used  to  give  an  analytical  expression 


which  relates  <AE>  i.  and  <AE>, 
all  d 


1 


Introduction 


The  collisional  vibrational  deactivation  efficiency,  p^,  has  been  used 

successfully  to  parameterize  thermal  unimolecular  reactions  involving  weak 

12  3 

colliders.  Tardy  and  Rabinovitch  '  and  Troe  have  developed  expressions 
relating  to  average  energy  transferred  quantities,  and  <AE>^j|, 

respectively.  These  quantities  are  defined  as: 

E  E 

<AE>^  =1  {E'-E)P(E',E)dE'/  f  P{E',E)dE' 

00  00 

<AE>^,l=  I  (E'-E)P(E',E)dE'/  |  P(E',E)dE' 

00  00 

also,  <AE>„  =  f  (E'-E)P(E',E)dE'/  f  P(E',E)dE' 

u  •'E 

where  P(e',E)  is  the  probability  that  a  molecule  with  internal  energy  E  will 
have  e'  after  collision;  additionally, 


where , 


-  <P>ap  <''E>p  -  <'»<io.n  <‘E>d  • 

00  00 

<p>  =  P(E',E)dE'/  I  P(E',E)dE' 

up  Je 

<'»do.p=  rP(E',E)dE'/fp(E'.E)dE' 


These  quantities  are  related  by  completeness,  P(E  ,E)dE  =1,  and  detailed 

balance,  (P  (E  '  ,  E) /P  (E,  E ' )  I //)^  exp  [- (E ' -E) /RT] )  .  Note  that  the  sign 

convention  makes  <AE)  ,,  =  -<AE>,  when  <p>.  =  1,  and  <AE>  ,,  =  <AE>  when 

all  a  down  all  u 

<p>do^n  ”  Although  the  P(e',E)  are  not  known  a  priori,  various  limiting 
model  types  have  been  used  for  comparative  purposes  (step  ladder,  Gaussian, 
exponential,  etc).^ 

2 

Some  confusion  has  resulted  because  <AE>,  and  <AE>  ,,  have  both  been 

d  all 

4-7  4 

used  in  the  literature  Recently,  Gilbert  has  shown  that  <CAE>^  is  the 

more  appropriate  quantity  to  use  in  parameterizing  p^ .  Oref  and  coworkers  in 
a  series  of  papers^  have  given  exhaustive  treatments  of  these  parameters  and 
their  interrelations,  for  the  case  especially  of  Bo  I tzmann- I i ke  transition 


2 


I  6 

probabilities.  For  a  given  P(E  ,E)  model,  Barker  has  also  shown  that 

<AE>^II  can  change  sign  with  temperature  even  though  <AE>^  does  not. 

In  the  present  work,  we  illustrate  dependence  of  the  ratio, 

7=<AE>^ I  I /<AE>^,  on  the  quantities  <AE^>  and  temperature  and  on  reactant 

properties  i.e.  internal  molecular  parameters  and  critical  energy  for 

reaction,  E  . 

o 

Ca I cu 1 ati ons 

Four  prototype  model  reactions  were  used;  decomposition  of  nitryl 
chloride  and  the  i somer i zati ons  of  methyl  isocyanide,  cyclopropane  and 
eye  I oheptatr i ene .  These  reactants  reflect  differences  in  vibrational 
frequency  patterns  and  molecular  complexity,  excitation  levels  and 
temperatures.  Table  1  exhibits  some  pertinent  quantities. 

For  simplicity  of  exposition,  a  step  ladder  model  of  the  coMisional 
transition  probabilities  with  step  size  AE  was  used  (AE  =  <AE>^  =  <AE>j) ,  and 
the  probabilities  were  computed  at  E=Eq' 

Pup  /Pdown  =  )exp(-AE/RT) 

0  0 

It  is  seen  that  the  general  definition  of  the  energy  ratio,  7=<AE>^ I  I /<AE>j, 
takes  the  form  of  eq.  1  for  a  step  ladder  model  and  is  a  function  of  step 
size,  temperature,  excitation  level  (critical  reaction  threshold)  and 
molecular  complexity;  the  last  two  dependencies  enter  through  the  density  of 

7sL  =(Pup-Pdown)  =  -  Pup/Pdown^/t^  *  Pup^Pdown^ 

=  -.AE^^E  +AE^^E  )®>'P(-*E/RT))  (1) 

0  0  o  o 

A  plot  of  7g|^  ^  (p^  +AE^^E  ^  ^’'P  ^  cange  of  temperatures  (200- 

°  -1 

16000  K)  and  step  sizes  (100-3200  cm  )  is  shown  in  Fig.  1  for  the  prototype 

systems.  As  required,  a  universal  plot  is  observed 

Equation  (1)  is  simplified  by  replacing  the  density  cf  states  for  the 

molecule  by  the  density  for  s  classical  oscillators  Define  p  =(p  /p  ,  )  , 

^  r  '"^up'  down^cl 


=  [1  .  j  e)cp[-AE/RT]  >  0; 


so,  75,  =  -(1-P^)/(1+P^)  • 

For  the  case  AE/E^  <<  1, 

Pr  =  (if  tl.(s-I).E/EJe-''^/«T  .  (2) 

o 

Thus,  7gl^  <  0  for  p^  >  1;  or  7g|^  >  0  for  p^  <  1  . 

These  limits  can  be  attained  at  high  and  low  temperatures,  respectively. 

In  the  low  temperature  limit,  the  Boltzmann  factor  dominates  so  that  7g|^  is  a 

function  mainly  of  AE  and  T.  In  this  limit,  increases  from  -1  with 

decrease  in  AE  and  with  increase  in  temperature.  In  the  high  temperature 

limit,  i.e.  exp (-AE/RT)'“1 ,  the  density  ratio  dominates  and  is  a  function 

of  AE  and  (s-l)/E^;  decreases  from  unity  with  decrease  in  AE. 

At  the  inversion  temperature,  Tj,  the  Boltzmann  and  density  factors 

compensate  one  another,  so  that  Pyp“PjQ^P  snd  p^=l .  Tj  is  sufficiently 

greater  than  typical  reaction  temperatures  so  that  it  is  only  approached  for 

reactants  having  an  unusually  small  E^  and/or  a  large  value  of  s.  For  the 

experimental  systems  under  examination,  this  temperature  might  only  be 

realized  for  eye  I oheptatr i ene  (Table  1).  Tj  is  the  temperature  for  which  E^ 

is  coincident  with  the  maximum  in  the  classical  energy  distribution.  It  is 

evident  from  eq .  (2)  that  Tj=E^/ [ (s-l)R] .  Thus,  to  this  approximation,  Tj  is 

independent  of  step  size  (AE) .  Indeed  for  the  prototype  systems,  only  a  weak 

dependence  on  AE  is  exhibited  in  Table  1.  Also  to  be  noted  from  Table  1  is 

that  a  larger  reactant  having  a  high  E^  shows  similar  values  (i.e.  density 

ratios  and  T.)  to  a  smaller  reactant  having  a  low  E  ,  eg.  C_H,_  and  CH„NC  . 

1  ^0^363 

When  +AE^^E  only  a  weak  function  of  molecular  complexity,  i.e.  for 

°  -1 

(s-l)AE<<E^,  (see  Table  1  for  the  case,  AE  <  400  cm  ),  7^^^  is  expected  to  be  a 
universal  function  of  AE/RT.  To  take  the  general  effect  of  molecular  complexity 


into  account,  an  efiective  temperature  can  be  used.  Tardy  and  Rabinovitch  have 
done  this  by  using  a  dimensionless  parameter,  <AE>/<E*>,  where, 


<E*>  =  J  (E-E^)p(E) 


-E/RT 


dE 


,(E)  dE 


lAe  now  define  an  effective  temperature,  as  the  temperature  at  which 
p^=exp(-AE/RT^);  so  that  =  T/ (1- (s-l)RT/E^)  =  T/(1-T/Tj)  =  T  Tj/(Tj-T); 
s  is  calculated  from  the  density  of  states  ratio  at  E^  i.e. 

— >-T 
e 

A  discontinuity  occurs  when  T  =  Tj ;  in  this  case  1/T^  =  0.  Figure  2  illus¬ 
trates  the  dependence  of  on  for  the  prototype  systems:  a  near-uni  versa 
curve  results.  A  shotgun  pattern  results  if  <AE>/RT  is  used  (not  shown). 


s=l+ln(/3^  +AE^^E  T<<Tj,  — ^>T,  and  for  T>>Tj,  ''”'1 


Earlier  work  in  the  literature  for  an  exponential  model  related  <AE> 


<ae^ii  -<ae>j 

and  <AE>^||  by  the  equation,  7^  =  =  ^AE^T  = 

a  d  a  c 

g 

The  Fj;  correction  was  added  in  later  work  and  was  stated  to  be  near  unity, 

g 

in  practice.  As  has  been  noted  ,  this  expression  does  not  take  into  account 

the  required  sign  change  for  7^  In  the  weak  collision  limit  (/?^  <0.03)  or 

<AE>^<<kT)  ,  7^j^p  -  -<AE>j/F^kT,  i.e.  is  linear  with  <AE>j .  For  the 

region  ( I  7^1^  I  <0 . 4)  ,  the  present  work  illustrates  that  7^|^  is  linear  with 

<AE>j  with  a  slope  of  ~  -1/(2RT^)  (Fig  2);  the  slope  depends  on  T^,  which  is 

usually  positive  Additionally,  if  eq  (3)  is  adopted  for  the  step  ladder 

case,  then  it  is  seen  that  F^  =  2[5 + (s-2)RT/E^]  or  F^>2,  always. 

Using  T^,  eq  1  giver.  7^j^=-tanh  (AE/2RT^)  .  Thus  we  ha .  e  determined  a 

relation  which  is  valid  for  all  step  sires  and  all  temperatures  and  which  is 

4 

substantially  easier  to  use  than  solving  an  1 ntegro-d i f f erent 1  a  I  equation 

From  the  above  relation,  <AE>  ,,  is  easilv  computed  from  <AE>  ,  and  T  Due 

'all  d  e 


to  the  transcendental  nature  of  the  furiction,  indirect,  methods  are  required 


to  obtain  <AE>^  from  I  I  and  (Regula  Falsi,  Newton-Raphson  or,  simply, 

by  iteration  on  Fig  2). 

In  future  work,  we  will  present  the  results  for  the  exponential  model 

Q 

which  can  be  directly  compared  to  the  expression  developed  by  Troe  ; 
different  analytical  expressions  for  will  also  be  evaluated. 

Cone  I  us i ons 

These  calculations  demonstrate  that  for  low  temperature  and  large  step 

size  I  I  ~  step  size  decreases,  or  temperature 

increases,  <AE>^ I  I  increases.  An  inversion  temperature  Tj  is  reached  when 

^down~^up'  condition  <AE>^ll=0.  With  increasing  temperature, 

Pup  >  P(jown'  ^^^^all  positive;  in  the  limit  <AE>^  I  I  =<AE>  ^ .  In 

general,  Tj  is  much  larger  than  normal  reaction  temperatures;  Tj  decreases 

with  molecular  complexity  and  increases  with  E^.  Thus,  for  large  molecules 

with  low  E  ,  the  difference  between  -<AE>  ,,  and  <AE>,  will  be  large, 
o  a  I  I  d 

The  functional  dependence  of  the  density  of  states  on  excitation  level 

and  number  of  oscillators  has  been  used  to  show  that  a  large  molecule  with  a 

high  E  behaves  in  a  manner  similar  to  a  small  molecule  with  low  E  ;  and  to 
o  o 

develop  a  universal  equation  relating  to  <AE>^  and  T^. 
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Table  1  Reactant  Parameter- 


FIGURE  CAPTIONS 


1.  Plot  of  vs  {p^  Ip^  )  for  nitryl  chloride,  methyl  iso- 

o  o 

cyanide,  cyclopropane  and  eye  I oheptatr i ene  for  step  sizes  of  100,  200, 
400,  800,  1600  and  3200  cm‘^,  each  at  250,  500,  1000,  2000,  4000  and 
8000  K. 

2.  Plot  of  7^1^  vs  <AE>^/RT^;  see  legend  of  Figure  1.  calculated  from 

expression  given  in  text. 
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